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Abstract 



This work addresses more to the technical rather than to the physical problem, how 
~f~, | to calculate analytically the form factor F(Q), the associated mean-square radius 

< r 2 >, and the distribution function <&(x, Q 2 ) for a given light-cone wave function 

^qq(x,k±) of the pion. They turn out to be functions of only one dimensionless 
pq | parameter, which is the ratio of the constituent quark mass and an effective Bohr 

momentum which measures the width of the wave function in momentum space. 

Both parameters are subject to change in the future, when the presently used so- 
^^ \ lution for the over simplified ||-model will be replaced by something better. Their 

relation to and agreement with experiment is discussed in detail. — The procedure 

can be generalized also to other hadrons. 
^ ^_ 

a: 
6-;, 

c-| , 1 Introduction and Motivation 



A quantitative measure of hadronic sizes is the mean-square radius. Its exper- 



X 

imental value for the pion (tt + ) is [1] J(r 2 ) = 0.67 ± 0.02 fm. One determines 
it by first measuring the electro-magnetic form factor F(Q 2 ) for sufficiently 
small values of the (Feynman-four-) momentum transfer Q 2 = — (p e — p e ') 2 , 
and then taking the derivative at sufficiently small Q 2 , i.e. 



{ r 2 ) = -6 dF ^ 



dQ 2 



Q 2 =0 



The electro-magnetic form factor can also be calculated. One of the most 
remarkable simplicities of the light-cone formalism [2] is that one can write 
down an exact expression. As was first shown by Drell and Yan [3], it is 
advantageous to choose a special coordinate frame to compute form factors 
and other current matrix elements at space-like photon momentum. In the 
Drell frame [4], the four- momentum transfer is —%cf = Q 2 = q 2 . The 
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space-like form factor for a hadron is just a sum of overlap integrals anal- 
ogous to the corresponding non-relativistic formula [3]. The general formula 
in [2] holds for any composite hadron and any initial or final spins S, but 
is particularly simple for a spin-zero hadron like a pion. It works with the 
wave functions \l/ n = ^ q q, ^ q q g , ■ ■ ., which are the Fock-space projections of 
the hadrons eigenstate. The total wave function for a meson is for example 
I* 'meson) = T,i(®qq(xi, k± t , \)\qq) + * 'qq g (xi, k ±i , \i)\qqg) + ■■■)■ The computa- 
tion of these projections is the goal of the light-cone approach to the bound- 
state problem in gauge theory [2], by solving Hlc\^) = M 2 |\P), with the 
eigenvalues M 2 being the invariant mass-squares of the physical mesons. 

It is advantageous to know Z 2 , the probability amplitude for finding the two- 
particle Fock-state \l/ u j in the pion total eigenstate \^ P i on ), see f.e. [2]. It can 
be obtained by computing the leptonic decay of the n + [4,2] as shown in Fig. 1, 

[ dxd 2 k ± r U AT j 16tt 3 , in [4]; 

Z2 i^T^ W = V^' W - = lVi<S?, present P) 



The factor N cn depends on the continuum normalization of the pion wave 
function, and n c = 3 is the number of colors. In Brodsky's convention, the 
empirical pion decay constant [5] is /„- ~ 130 MeV/-\/2 rs 93 MeV. The rela- 
tion involves only the L z — S z — component of the general ud wave function, 
where i/)(x, k±) = ^ u ^(x, k±; ]{)/Z 2 is the (normalized) probability amplitude 
for finding the quarks with anti-parallel helicities, particularly for finding the 
up-quark with longitudinal momentum fraction x and transversal momentum 
k±, and the down antiquark with 1 — x and —k±. Their respective charges are 
ei and e 2 , respectively, with e± + e 2 — 1. 

The analytical calculation of Z 2 is the second aim of this work. As shown 
below, Z 2 is not only finite but even large, see for example [6,7]. The third 
aim of this work is to calculate the pions distribution function 

Q 2 

®ud(x, Q 2 )=k J dk 2 ± Vjfc, k ± ; U), (3) 

o 



where the transversal momenta are integrated up to some momentum scale 
Q 2 [2,4]. The distribution function continues to play an important role and we 
cite only a small fraction of the available literature [4,8-15]. 

It appears as if this paper is loaded with formalism, but an attempt was made 
to proceed pedagogically. In section 2, we expand on general considerations 
and give the explicit formulas for the purpose of definition and notation. In 
section 3, the same results as in a precursor to this work [18] are derived in 
a simpler way. The non-interested reader may proceed to section 4, where 
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Fig. 1. Matrix element of the charged 
axial-vector current controlling the de- 
cay 7T — ► ev. 



Fig. 2. The reduced wave function (f(p) 
for the pion is plotted versus p in an ar- 
bitrary normalization. The filled circles 
indicate the numerical results, the open 
circles the fit to Eq.(8). 



the necessary integrations of Eqs.(l-3) are carried out explicitly quoting only 
the definitions and the results. All intermediate steps are left out, albeit they 
contain the lion's share of this work. The emphasis of the present work is on 
the analytical evaluation of expressions taken from the literature, but some 
physical aspects are discussed in section 5. In section 6, explicit expressions of 
the distribution function are discussed, and finally, in section 7, the necessarily 
compact presentation of this work is summarized. 



2 General considerations 



It is known empirically, that the form factor at low Q 2 has essentially mono- 
pole structure [1], and that the mean-square radius is essentially all the infor- 
mation there is. In the sequel, we restrict considerations to the contribution 
of the 2-particle Fock-state to the form factor of the pion, i. e. to 

F2{Q) = dx d 2 ki_(eiip{x,k±_ + (1 - x)q±) + e 2 ip(x,k± - xq±)J ip(x,k±).(4:) 



Since ip(x,k±) is normalized, F 2 (0) = 1. The associated mean-square radius, 
the 'size', is denned by 



R\ = (r 2 ) 2 = -6 



dF 2 {Q 2 



dQ 2 



(5) 



The three quantities R 2 , Z 2 and $ U( j(a;, Q 2 ) are of physical interest, since they 
can be measured [1,5,16], to some extent, see also section 5. For a given light- 



cone wave function ip(x, k±) their calculation is straight-forward. But here 
is a problem: The bound-state equation for non-perturbative QCD, which 
should define this function, continues to be a challenge [2]. Recently, how- 
ever, some progress was made with the oversimplified fj-model [20] for the 
pion and other pseudo-scalar mesons. It produces a numerical ip(x,k±). But 
the three-dimensional numerical integration required in Eq.(4) and its subse- 
quent derivation with respect to Q 2 is cumbersome and may be numerically 
inaccurate. It might be (numerically) more reliable to suitably parametrize 
ip(x,k±) and to perform the quadratures invoked by Eq.(4) analytically. In 
fact, a certain form of the parametrization is imposed by the structure of the 
bound-state equation in general, see however also the precautions mentioned 
in Section 7. 

Quite in general, one is able to write down an integral equation in the three 
variables x and k± for the wavefunction ip(x, k±), see f.e. [2,20]. The solution 
of such an equation is numerically non-trivial, among other reasons, because 
the longitudinal momentum fractions are limited to < x < 1. It is therefore 
advantageous to substitute the integration variable x by an other integration 
variable k z , — oo < k z < oo, which has the same range than either of the 
two transversal momenta k±. For equal quark masses m u = m^ = m, the 
substitution by the generalized Sawicki transform is simple, see for example 
[21], i.e. 

*(fc«) = ^|l+ 1 = 1 «=* *2 = ("» 2 + ^fe^-( 6 ) 



Formally, the three integration variables k z and k± look like a conventional 
3-vector p= (k z ,k±). The substitution and the associated Jacobian seems to 
destroy the hermitian property of the kernel, but only superficially, since it 
can be restored by substituting 




v/l + p 2 /m? 
ijj(x, k±) = ip(k z , k±) (7) 

Jx(l — x) 



Mathematically, the so obtained integral equation for the reduced wave func- 
tion <p(p) is identical with the original integral equation for the light-cone 
wave function ip(x, k±). It looks like an integral equation in usual momentum 
space (p), but continues to be a relativistically correct front-form equation. 
The reduced wave function ip(p) describes a bound state and must decrease 
for p — > oo with a certain scale p a , either exponentially or like a power. As dis- 
cussed in [18] a power law is more likely. In the present context it can suitably 



be parametrized by 

V(P) = n ~ 2/ 2f (8) 

The analogue to the Bohr momentum p a can be fitted to the numerical solu- 
tion. An example is given in Fig. 2. More details can be found in [2,19,20]. 

To be even more specific and concrete, the integral equation in the jj-model 
is taken from [19] as an illustration without proof, i.e. 



M 2 ip(p) = Am 2 + 4p 2 ip(p) 

a f d 3 p' ( Am 2 2/i 2 



3 2tt 2 J m \(p — p') 2 jj? + {p — p 



r\2 



tptf). (9) 



It looks very simple, indeed, and is rotationally invariant. All well-known diffi- 
culties of the front form with rotational invariance seem to be absorbed by the 
factor \/x(l — x) in the mapping Eq.(7). In [19] it has been shown how to solve 
the equation numerically for spherical symmetry tp{p) = <p(p)- Some further 
results are compiled in Table 1, below. Here and below, masses and momenta 
are expressed in units of u — 350 MeV, except when noted otherwise. All 
parameter sets produce a lowest mass eigenvalue of M 2 = (139.57018 MeV) 2 , 
the mass (squared) of the n + . Both the M 2 and the two first exited states 
are stationary with respect to /i. The (unphysical) regularization parameter \i 
determines itself from the solution ('renormalization'), for details see [19]. 

Next, before proceeding with the computation of the form factor, a number of 
notational definitions are introduced, in terms of which the final results turn 
out to be simple. Once one has ip(p) in a parametrized form like Eq.(8), one 
can transform back to the variables x and k± and define 

Z(x,k ± ) = l + £ = l + ^^f^ + — 2 )\ y (10) 

Pa Pa X \ l - X ) ^Pa X \ l ~ X ) 

as well as X(x) = Z(x,k± = 0). The combination 1 +p 2 /m 2 is trivially 
obtained from this by putting p a = m. The dimensionless parameters s and t, 

m 9 m 2 — p 2 

s = — , and t 2 = r^, (11) 

Pa rn 2 



will govern the results below. The form factor in Eq.(4) will be calculated in 
the form 



1 oo 



F 2 {q[)= J dx J dkl\^j{x,k ± )\ 2 f{x,k ± ;q±), (12) 





2tt 



, ei^(x,fcj. + (l - x)gL) + e 2 ^(a:, k ± - xql) 

with /(x, fej.; a.) = / d<f> — — — =- • (13) 

J Q 2^(x,k jL ) 

The function f(x, k±;q±) contains all the difficulty in the problem, particularly 
the integration over the angle between k± and q±. It depends only on the 
absolute value of the momentum transfer. If one restricts to calculate only the 
derivative with respect to q±, as in 





1 oo 

(r 2 } 2 = [ dx f dkl\ip(. 




v,k±)\ 2 g(x,k±), 




(14) 


with 


g(x,k ± ) = -Q f(x, k±;q±) 


n = - 6 -r^f(x,k±;q±) 


1 
o,=0 


(15) 



one deals with a function g(x, k±) which is independent of q±. 

For calculating the probability amplitude according to Eq.(2) with a normal- 
ized wave function as in Eq.(7), and the mean-square radius according to 
Eq.(14), one needs to evaluate three integrals, namely 



Air r 

—p 2 a N{s)M 2 = j dxd 2 k ± ^ 2 (x, k ± ), ( = 1 ) 

2-Kp 2 a P(s)N = dxd 2 k ± if) (x,k±), 

ttS(s)N 2 = f dxd 2 k ± ^ 2 {x, k L )g{x, k±). ( = (r 2 ) ) (16) 



They are cast into three dimensionless functions N(s), P(s), and S(s), which 
depend only on the dimensionless parameter s as defined in Eq.(ll). Once 
they are known, the calculation of the probability amplitude and the root- 
mean-square radius is easy. In particular, with 



_ U 2ir ^N(s) _ 0.5665 JW) 
2 Pa 3 P(s) m S P{s) ' [ ' 



R _V3 ^S(s) V3 

K 2-~ 1 , tie — ~ , U 

2 Pa yjN(s) 2 Pa 

one can study them as dimensionless functions of s for a fixed value of m. 



3 An approximate treatment 



As will be seen below, the factor Jl + p 2 /m 2 in Eq.(7) poses certain problems 
in the evaluation of the form factor and other integrals. In order to proceed 
pedagogically, this factor is first replaced here by unity, by using 



ip(x,k±) = 



N 



x(l 



X) 



'1+P 2 /Pl) 



(19) 



with the two adjustable parameters k and p a . This is justified, of course, only 
if the Bohr momentum obeys p 2 -C m 2 , i.e. in the non-relativistic limit. But 
here p a /m ~ 1.2: The quarks move highly relativistically! 

We disregard this objection, and proceed. For to compute the mean-square 
radius, one needs to evaluate the function / as defined in Eq.(12). For the 
parametrization of Eq.(19) one has thus 



f(x,k ± ;q ± ) 



2l ZK 



2~ 



+ fz« 





2tt 



+ 



m 2 {x - \) 2 , {k±~(l-x)q L 
Vl x(l - x) 



i | m 2 (x - \) 2 | 



Ap 2 x{l — x) 



k± — xq± 



p 2 x(l — x) Ap 2 a x{l — x) 



It has two contributions, one from the quark (~ e\) and one from the anti- 
quark (~ e<i)] they differ from each other by exchanging x with 1 — x. Since 
one is interested in the ^-integration, one rewrites the equation as 



71 

f(x,kr,q±) = Z K J l 



( : 'i 



{pi + C\ COS ( 



^ + 



e 2 



(b 2 + c 2 cos 4>) K 



(20) 



where the coefficient functions b and c for the quark are given by 



b x = Z 
6i 



1 — x ql 

x Ap 2 a ' 
1 1-x 



4p 2 



x 



<h 



CiCi 



1 k± q_L 

X Pa 2>Pa ' 

1 1-x 



2 Pa X 



(Z-X). 



(21) 



f22) 



The functions Z(x, k±) and X(x) had been defined in Eq.(10). The coefficients 
for the anti-quark, i.e. b 2 , c 2 , b 2 and c 2 , are obtained from those for the quark 
by exchanging x with 1 — x. 



However, one need not know explicitly f(q±)- As a big advantage of the present 
aim to calculate (r 2 ), one needs the form factor only for very small qj_. Since 
C\ and Ci vanish for qj_ — > 0, one can expand Eq.(20) with c before integration. 
The contribution from the quark becomes then 



Z K r 
fi(x,k±;q±) = ei— J 



d k(k+1) c\ 2 
1 — K— COS0H 777 COS « 

t>i 2 of 



+ 



.(23) 



Note that the term of first order in C\jh\ vanishes upon integration, and that 
the term of second order is integrated trivially. Taking the derivative with 
respect to qj_ according to Eq.(15), one is left with 



g 1 (x,k±) =<oe 1 ix 



Z K 
~bf 



b\ k(k + 1) c\b\ 



h 



hi 



After inserting the derivatives according to Eq.(22) one gets 



gi(x,k± 



3ttk 
2^" 



K 

~z 



(«+l) 



X 

~z* 



x 



ei- 



x 



(24) 



The total g — g\ + g-i is thus proportional to 



1 — x x [(ei + e 2 )x 2 + ei(l — 2x)] 

ei he 2 - 



,r 



'1 — x 



x(l — x) 



X 



x(l — x) 



The great simplification occurs since e x + e 2 = 1, and since both X and Z 
are symmetric under the exchange of x and 1 — x, such that 1 — 2x vanishes 
upon integration. Finally, the (^-function of Eq.(15) for the semi-relativistic 
wave function, Eq.(19), becomes for arbitrary k 



g(x, k 



■± 



3ttk 
2^ 



K 

~z 



X 



x 



x(l 



X 



(25) 



One should emphasize that g(x, k±) was obtained without explicitly evaluating 
Eq.(20) and that the results for k — 2 and k = 3/2 agree with the expressions 
in [18], where the limit q± — > had been taken after integration. 

Having g from Eq.(25) and the wave function from Eq.(19), one can calculate 
the three integrals as defined in Eq.(16), particularly 



3 1 



N(s) = —j^Jdxdkl^(x,k ± ), P(s) 



2plN 



dxdk]_ ip(x, k± 



1 1 f 

S ^ = ~Tf2 J dxdk l V> 0&, k±) g(x,k±). (26) 



Before evaluating them, it is advantageous to change the integration variables 

kl 



x and k\ to z and u, respectively, with 



k 2 

z = 2x-l, and u=^r. (27) 

m 2 



Z(x, k±) and X(x) as defined in Eq.(lO) are then rewritten as 

Z(,,») = ,» " + (; = y + "' , X(,) = Al±fcW. (28) 



Evaluating Eq.(26) for k — 2 gives the normalization function 



Ms) = — / — / du - ^ t 29 

V ; <fi I o. I [l + (z 2 - l)t 2 + uf y ! 



-1 

-8 - 10s 2 + 3s 4 + 36(s)s 2 (8 - 4s 2 + s 4 
8(s 2 -l) 3 



(30) 



Here and below the abbreviation b(s) = arctan(V s 2 — 1)/ \Js 2 — 1 is used. The 
major labor of the present work is hidden in the evaluation of integrals such as 
Eq.(29) and to express them as explicit functions of s. In order not to load the 
paper with straightforward formalism, the arithmetics are suppressed here, 
as a rule, the more as they can be produced also by Mathematica. Care was 
taken that no typos do occur. -- For the probability function one obtains 

„, , 1 fdz 7, [l-^ 2 ] 1 tt1-3s 2 + 2s 3 . . 
P( s ) = — — du - - y = — n ^—, 31 



and the size function finally is 



12 r dz 



1 , oo 



s 8 J 2 

-1 



l2 



S(s) = - - du[l + zY 1-z 2 ' , ' ■ l6 J (32) 



3 [1 + {z 2 - l)t 2 - 2«] 



[1 + (z 2 - l)t 2 + u] 
-136 + 72s 2 - 56s 4 + 15s 6 + 36(s)(32 - 16s 2 + 36s 4 - 22s 6 + 5s 6 



40(s 2 - l) 4 

These functions are plotted below as function of s in the region of interest 
s ~ 1. For k — | they are given in [18]. 



4 The full treatment 



Having been so explicit in the previous section, we can be more compact when 



restoring the factor Jl + p 2 /m 2 in the wave function 

~, (rM - ; £ y/l+P 2 M 2 _ M YHx,k ± ) 

* { ' ] ~ fi^x) (i + pVpD 2 " v^I^) ^ k± ) ■ [66) 



The origin of the coefficient function Y(x, k±) is obvious, i.e. 

Y{x,k A _)^l + — 2 = A n t^ kl v Y ^ U ) = TT ± 4- (34) 
m 2 4m 2 x(l-x) [1 - z 2 ] 



Correspondingly one has U(x) = F(a:, /cj_ = 0) and C/(z) = Y(z,u = 0). The 
coefficient function Z was defined in Eq.(10). Note that the full relativistic 
wave function interpolates to some extent between k — 2 and n — § in the 
previous section: Substituting m by p a gives the same wave function as for 
k — |, while putting m — > oo gives the wave function for k — 2. 

For to compute the form factor, one needs to evaluate / as defined in Eq.(12). 
The analogue of Eq.(20) is 



Z 2 f 

Y 2 J 



1 1 

2 



(Si + d COS 0) 2 (Sa + C 2 cos , 

e i — — : — + e 2^ — ; 7T2 



.(35) 



{pi + d COS 0) J (6 2 + c 2 cos 

The coefficient functions b{ and q had been given above, and 



B _ Y | l ~ x gi ^i = -i — — (36) 

x 4m 2 ' a; m 2m ' 

j 1 O" . II — x 

Bi = ^~ 2 -, C X C X = — (Y-U). (37) 

4 m z x 2 vn A x 

Again, the anti-quark coefficients are obtained by exchanging x with 1 — x. 

The analogue of Eq.(23) is now much more complicated but still straightfor- 
ward. In full analogy with Eq.(25) one gets finally 



g(x, k±) = -j- 

Pa 



2 3X 1 / 7 X U 



Z Z 2 s 2 \8Y ZY 8Y 2 



„2 



(3* 



x(l — x) 



Note that this equation is consistent with the previous results: Omitting the 
terms in the round bracket (thus m — > oo) agrees with Eq.(25) for k — 2, 



10 



while formally putting s = 1 and Y = Z reproduces it for k, 



2' 



Having g and the wave function, one can evaluate the three integrals as defined 
in Eq.(16). The normalization function becomes 



3 fdz 7 . \l-z 2 } 2 \l + u] 



N ( s ) = ~, /-*-/<*« 



s 6 _7 2 V [1 + (z 2 - l)t 2 + u 

4 - 4s 2 + 3s 4 + 36(s)s 4 (-2 + s 2 ) 
8s 2 (s 2 - l) 2 



i 



(39) 



and the probability function is 



, v _ + fdz J du [1-z 2 ] [1 + ufi s 2 + b(s)(s + s 2 -s 3 -2s 4 ) 
(SJ "i 2 {[l + {z 2 -l)t 2 + uf" 2(1 -O -^ 



Its integration is non-trivial, see App. A. The size function 5'(s) is now much 
more complicated than in the previous section, but after elementary manipu- 
lations one gets 



12 fdz 7 . r ,9 r 9i2 / 2fl + «l 

1 - z ' 



s{s) = ^J Y J du [1 + z 

-1 o 



2 



[1 + (,2 2 - l)t 2 + U 



5 



3[l+ M ][l + (^ 2 -l)t 2 ] , 7 



+ — r4 ^ ^ + 



[l + (^ 2 -l)£ 2 + «] 6 8[1 + (^-1)£ 2 + m] 4 

[1 + (z 2 - l)t 2 ] 1 



[l + (^ 2 -l)£ 2 + «] 5 8[l + «] [1 + (z 2 - l)t 2 + u]\ 



(41) 



_ 34 + 37s 2 - 41s 4 + 15s 6 + 36(s)(-8 - 16s 2 + 21s 4 - 17s 6 + 5s 8 ) 
~~ 40s 2 (s 2 - l) 3 ' 

This completes part of the aim of the present work, namely to calculate ex- 
plicitly the mean-square radius (r 2 ) and the probability amplitude Z 2 , see 
Eqs.(17) and (18). 



5 Discussion 



In a non-relativistic system one calculates the mean-square radius compara- 
tively cheaply by Fourier transforming a momentum space wave function to 
configuration space, and taking the expectation value of r 2 . Taking for sim- 
plicity the function <f(p) in Eq.(8) as such a 'wave function' with all due pre- 
caution, precisely this was intended in [19], with the resulting non-relativistic 
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Fig. 3. The discrepancy d = R e jRi 
is plotted versus s by the solid line. 
It is compared to the semi-relativistic 
case by the dashed (k = 2) and the 
dashed-dotted line (k = 3/2). 




Fig. 4. The rms-radius of the pion in 
units of 0.67 fm is plotted versus s by 
the solid line. It is compared with the 
approximate cases by the dashed (k = 2) 
and the dashed-dotted line (k = 3/2). 



estimate R e = \^3/2p a given in Eq.(18). It hurts to admit, that a missing fac- 
tor 2 in [19] was unraveled only in the course of the present work, particularly 
in [18]. But now, we are in a better shape: Not only are we able to calculate the 
mean-square radius in the almost same way as an experimentalist performs 
the measurement, but we can calculate it analytically, irrespective of whether 
one deals with a non-relativistic system or not. How does the discrepancy 
behave? 



The discrepancy d = R e /R.2 = yN(s)/S(s) is plotted versus s = m/p a in 
Fig. 3. One observes that the curve for the approximate wave function of 
Eq.(19) for k = 2 approaches the limiting value d = 1 rather quickly from 
below for growing s. Note that there is a long way to go until s ~ 137, the 
value for a Bohr atom with Bohr momentum p a ~ ma. This can also be seen 
from the analytic formulas, Eqs.(29) and (32). The case for k = 3/2 has also 
a limit, however uninteresting in this context. Fig. 3 displays as well that 
the correct wave function from Eq.(33) approaches the limit from above. The 
discrepancy exceeds rarely some few percent for values of s down to s ~ 0.3. 
That the non-relativistic estimate is so accurate was a surprise. Ultimately, 
for p a — > oo thus s — > 0, for a point like system, the discrepancy and thus the 
size drops to zero even faster than the non-relativistic estimate. This can be 
seen also in Fig. 4, where the size is plotted for the various cases considered 
in this work. 



Table 1 summarizes the results for solving Eq.(9) numerically, in the manner 
described in [19]. Set #1 quotes the results of [19]. The entry R e = 0.33 fm for 
the non-relativistic estimate reflects the factor 2 mentioned. The calculation 
of the exact root-mean-square radius i?2 = 0.30 fm was one of the motivations 
to get this work started. The discrepancy d = 0.92 is amazingly small. Two 
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Table 1 
Compilation of results. 
















set a m 
MeV 


Pa 

MeV/c 


Pa/m 


s 


d 


z 2 


Re 

fm 


R 2 

fm 



#1 0.690 406 515 1.268 0.789 0.919 0.392 0.33 0.30 

#2 0.682 301 378 1.255 0.797 0.919 0.524 0.45 0.42 

#3 0.672 196 234 1.192 0.839 0.920 0.774 0.73 0.67 

additional calculations have been performed, one with the constraint R 2 ~- 



{f 2 )exp — 0.67 fm in order to pin down the quark mass m, and the other 
with R 2 = 0.42 fm, referred to as set #3 and #2, respectively. The quark 
mass m determined from the constraint on the rms is almost linear to 1/R 2 . 
But the width of the reduced wave function changes almost in proportion, 
such that Pa/in ~ 1.2 (or s ~ 0.8) and thus d is almost independent of the 
case. No explanation for this numerical observation can be offered. 

Table 1 includes also the probability amplitude Z 2 , as calculated with Eq.(17). 
Early estimates of this quantity from Brodsky, Huang and Lepage [6,7] imply 
that the pion is roughly 1/4 of the time in the gg-Fock state, but it is the first 
time, that such a number is calculated from a wave function. The numbers in 
the table imply 1/2 to 1/5 for this quantity, depending on the case, in rough 
agreement with [6,7]. 

As a matter of fact, Brodsky, Huang and Lepage emphasize also that 2-particle 
Fock-state is smaller or more compact than the pion according to Eq.(l). The 
estimate they give in [7] (on pg. 20) is mutatis mutandis R 2 — Z 2 ■ (0.7 fm), 
where Z\ (not to be confused with the Z 2 of Eq.(17)) is the probability to 
find a bare quark in a dressed quark. It is certainly less than 1, but neither 
measured nor computed, thus far. One concludes that a fit to the electro- 
magnetic size as in set #3 is the least favorable of the three sets in Table 1. 
Set #1 corresponds to a 'hard core radius' of about R 2 ~ 0.3 fm, corresponding 
to a Z 2 ~ |. Some people favor a hard core radius of 0.4 fm, which was the 
motivation for calculating with set #2. 



6 The distribution function 



In a recent experiment, Ashery [16,17] has provided the first direct measure- 
ment of the pion light-cone wave function (squared): A high energy pion dis- 
sociates diffractivly on a heavy nuclear target. In a coherent process the quark 
and the antiquark break apart and hadronize into two jets. Their momentum 
distribution carries information on the quarks momentum distribution in the 
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pion. By and large, Ashery's results agree with the x(l — ^-distribution pre- 
dicted long ago by Lepage and Brodsky [4] and by Efremov and Radyushkin 
[8], by solving a perturbative QCD evolution equation in the limit of large 
momentum transfer Q 2 — »■ oo. Ashery has fitted his data by a linear superpo- 
sition with the x(l — x)(2x — l) 2 -distribution of Chernyak and Zhitnitsky [12], 
which is believed to be reliable in the opposite limit Q 2 — + 0. 



Although Ashery's experiment cannot be analyzed in terms of the distribution 
function alone, it is interesting to ask what the present work predicts for 
it. - By convenience, the definition in Eq.(3) is split up into a normalization 



factor Nq, and the reduced distribution function 
N^(f)(x,k 2 ), with 



(x,k 2 ) , i.e. $ g q(x,k 2 ) 



Ns, = AnMZ 2 ^ = f n 



4tt p 2 a R 



m^ P(s] m 2 V n, 

(x, k 2 ) = 2m x( 1 — x) / dk 2 

o 




[kl + a 2 ] 2 ' 
a 2 =m 2 + 4x(l — x) yp 2 a — m 2 ) . 



(42) 

(43) 
(44) 



Note that a 2 



a (X, 



actually is a function of x and that /„- is the carrier of 



the dimension. The factor of 2m is inserted into the definition of <p to have a 
dimensionless function. In Eq.(42), M and Z 2 have been substituted according 
to Eqs.(16) and (17), respectively. Brodsky has emphasized repeatedly [2] that 
one could normalize the wave function also by the weak pion decay amplitude, 
Eq.(2). It is gratifying to see that the normalization function N(s) cancels in 
Eq.(42) even without taking special care of that. Evaluating the integral in 
Eq.(43), one gets for p a > m 



(f>(x, k 2 ) =x(l — x) 



2m z 



+ 



2m 



and for p a < m 



yc 



m z 



2m 



arctgi 



[k 2 + m 2 } 
n + a 2 } 
I a 2 + k 2 



a 2 — m 2 



— arctgi 



m z 



(45) 



(x, k 2 ) 



x(l — x) 




k 2 + m? — \Jm 2 — a 2 



k 2 + m? + \Jm 2 — a 2 
m + y/m 2 — a 2 \ 
- y/m 2 — a 2 J 



(46) 
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Fig. 5. The Fock-state probability am- Fig. 6. The reduced asymptotic dis- 

plitude Z2 is plotted versus s by the solid tribution function (j) as ( x ) as given by 

line. It is compared to the approximate Eq.(47). The full line refers to p a = 3m, 

cases by the dashed (k = 2) and the the dash-dotted line to p a = rn, and the 

dashed-dotted line (« = 3/2). dashed line to p a = m/3. 

The quadratures are elementary: One checks them by differentiating the latter 
two equations with respect to k\. 

In the limit k t — > 00, one gets from this the asymptotic distribution function 



,(x) =x(l — x) • < 



2m 2 



2m 2 



+ 



+ 



2m 



y/c 



m^ 



=arctgi 



\fm? 



m 



2m , , m 2 
=ArthW — 



m* 



for p a > m, 

for p a = m, (47) 

for p a < m. 



One notes first that the Lepage-Brodsky-Efremov-Raduyshkin limit (LeBER) 
is reproduced by the overall factor x(l — x). For p a = m, we reproduce it 
even exactly. One should however not be surprised by deviations from this 
limit, since hadronic scales like m or p a had been omitted consistently in the 
derivation [4,8], as producing terms of order unity. On the other hand, such 
terms can generate significant deviations and sometimes can change even the 
qualitative behaviour, as will be discussed next. 

In the familiar bound-state problems, the constituents have a momentum 
whose mean (~ p a ) is significantly smaller than their mass. As shown in Fig. 6, 
such non-relativistic systems with p\ <C m 2 have an asymptotic distribution 
which is (much) narrower than the LeBER distribution x(l — x). In the ex- 
treme non-relativistic limit (as for example in a /ie-atom with equal masses), 
the distribution function will be peaked very sharply at x = \. The distri- 
bution functions will be plotted here in units of the peak value. They are 
functions only of s = m/p a , and not of m and p a separately. The reason is 
that the dimension of $ is carried by the pion decay constant, as mentioned. 
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Fig. 7. The distribution function 
(f)(x,kt) with p a = m/3 is plotted ver- 
sus x for k t = 32.6, 131, 262, 588 MeV 
by the four dashed lines, respectively. 
The solid line gives the asymptotic dis- 
tribution function (p as {x). All functions 
are plotted in units of the peak value 
M (i) = 5.435. 



Fig. 8. 



(x, ltf) with p a 



The distribution function 
3m is plotted versus 
x for k t = 0.647, 1.18, 2.35, 5.29 GeV by 
the four dashed lines, respectively. The 
solid line gives the asymptotic distribu- 
tion function (f) as {x). All functions are 
plotted in units of 4> a s(h) = 0.1892. 



We are less familiar with bound-state systems in which the constituents have 
a mean-momentum larger than the mass. Such systems are hypothetical. Let 
us refer to them as relativistic bound-state systems, characterized by p 2 a > m 2 . 
As displayed in Fig. 6, a (highly) relativistic system has a distribution function 
(much) broader than the LeBER limit. In particular, it is natter around x = |. 
In the limit p a — > oo, it becomes completely flat, which reminds one to the 
point-like Schwinger/t'Hooft bosons in 1-space and 1-time dimension [2]. 

Similar considerations hold true also for finite values of k t . In Figs. 7 and 8 
the distribution function <fi(x, k 2 ) is plotted for a non-relativistic and a rela- 
tivistic bound-state system, respectively, for different values of k t . One should 
emphasize that the relativistic case develops 'ears' for sufficiently small k t , 
which remind to the Chernyak-Zhitnitsky distribution. These are absent in 
the non-relativistic case. 



7 Summary and conclusion 



The essential progress of this work is the consequent application of the light- 
cone wave function ip(x,kj_) in the parametrization of Eq.(7), i.e. in terms 
of a reduced wave function (p(x,k±). The reduced wave function in turn is 
parametrized in Eq.(8) like a Coulomb wave function in momentum space. 
The only adjustable parameter is the analogue of the Bohr momentum p a . It is 
adjusted to a numerical solution of a mock-up, the Tj-model. Both steps com- 
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bined allow to calculate many observables, particularly the mean-square radius 
(r 2 ) 2 , the probability amplitude Z 2 , and the distribution function $ q q(x, k%) 
explicitly and analytically as functions of m and p a . The calculation of Z 2 is 
novel, particularly for QCD. 

The weak point of the present work is its model dependence. The tj-model 
picks out one particular aspect, namely the strong attraction of the spin-spin 
interaction in the singlet channel. While such description is perhaps viable as 
indicated by several other phenomenological models, the above discussion on 
the model results might be misleading. As listed in Table 1, the effective quark 
mass in Eq.(9) is much larger than the bare quark mass in the original QCD- 
Lagrangian. Thus, the solution of reduced wave function in Eq.(9) should be 
regarded as the bound state of not bare but dressed quark and antiquark. The 
effective 2-particle bound state is then no longer a compact object. In fact, 
the form factor obtained only with the subset (||) component of light-cone 
helicity amplitude cannot be trusted. In multiplying the two light-cone wave 
functions to compute the form factor, the || and J, J, components contribute 
perhaps as significantly as the || and || components. There are also other 
aspects which must be improved in the future. 

One should emphasize however, that the | j-model works only with the param- 
eters appearing in the QCD-Lagrangian, the quark masses and the coupling 
constant. This is kind of the minimum requirement for any effective theory, 
inspired by QCD, and an improved solution cannot have less information than 
the momentum spread. The model applied serves at least the purpose to relate 
p a to m and a through the (numerical) solution. 
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A On the evaluation of the integrals 



The major labor of this work is in the evaluation of the integrals. As a rule, they 
are straightforward, with one exception: Evaluating Eq.(40) as it stands, even 
Mathematica drives crazy. Therefore some subtle aspects should be mentioned. 
From the definition 

, , fdz 7 [l-z 2 ] fl + nl* , k N 
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follows the particular value -P(O) = 4/3. Introducing the new integration 
variable x = \/l — z 2 and integrating by parts over u, two terms arise in 
P(t) = Pi(t) + P 2 (t) with 

„ . . r dxx 3 1 „ . . f dxx 2 Arth.(xt) ., , 

™ = / vt=fi-w and ™ = I vT^ t ■ (A - 2) 

v v 



Note that -Pi(O) = ^(O) = 2/3. The first of them is elementary, i.e. 

. F\ Aril, I l/v/'l - r-i\ /,/'! - /- ! 

Pl(t) 



^T^2 Arth ((tVT^) jJY—p 

t 2 t 2 VT^¥ 



(A.3) 



while P2{t) has resisted all attempts to integrate it up. Only when observing 
the identity d(tP 2 (t))/dt = Pi(t) it was possible to proceed with 

P 2 (t) = - t JdyP 1 (y). (A.4) 



The (vanishing) integration constant can be determined at the end from re- 
quiring P 2 {0) = 2/3. One can insert either the power expansion for 

r, , , 2 8 t 2 16 t 4 r , , , A rN 

PlW = 3 + 5F3 + T^ + "" fOT |(|<1 ' (A ' 5) 



or one can integrate it up in terms of the useful function 



K(x) = - I — = 1 + t -f^B^ 2 " (A.6) 

xJ sin a; ~ L {2n + 1)1 

m^ IT 

where the B n are Bernoulli's numbers. One gets to this form from Eq.(A.4) 
by a final variable transform y = sin a;, i.e. 



arcsint 

1 1 f dxx 



W) = ¥ + l I 



sin 3 :r 



1 arcsint y/1 — t 2 arcsint . . . , 

= — — + F n (arcsint). (A.8) 

It 2 It t 2 2t y ' y ' 

In the approximate Eq.(40), F u { arcsint) was approximated by 1. 
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